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Description 

[0001] Thi3 invention relates to communication systems and particularly to Ihose employing orthogonal frequency 
division multiplexing (OFDM) modulation methods, 
s [0002] OFDM is a technique whereby information symbols are communicated from a transmitter to a receiver over 
a plurality of sub-carriers. 

[0003] Carrier and clock frequency offsets can produce large degradations of the Bit Error Rate (BER) performance 
of Orthogonal Frequency Division Multiplex (OFDM) systems . See for example Proc. IEEE 83(6) 982-996 June 1 995 
and Proc. COST254 1997, Indeed they not only produce extra-noise due to Inter Carrier Interference (ICI) but also a 
10 parasitic rotation of the symbols which also increases the BE R. 

[0004] Herein, carrier frequency offset means any difference in frequency between the carrier Irequency generators 
In the transmitting and receiving circuitry and clock frequency offset means any difference between the sampling rates 
of the transmitting and receiving circuitry. 

[0005] There already exist several methods for estimating and compensating a carrier frequency offset. See for 
w example, US-A-5450456 (Mueller) and Elec. Lett. 33(2) 113-114 Jan. 1997. However, using either of these known meth- 
ods, the residual carrier frequency offset can rise up to a few percent of the sub-carrier spacing and there still remains 
a degradation. As regards clock frequency offset, the degradation depends cm the number of sub-carriers. II this number 
is large, this may prevent the designer from relaxing constraints on the receiver's clock generation mechanism. 
[0006] To Illustrate the problem, consider a typical OFDM transmission system such as Digital Audio Broadcasting 
20 (DAB) mode 1. In fig. J and fig. 2, the degradation due to ICI is plotted. After correction with the classical algorithms, 
the carrier frequency offset can reach 3% of the sub-carrier spacing (in fig. 2, this corresponds to f= 3% where f is the 
carrier frequency offset expressed in number of sub-carrier spacings, N is the number of sub-carriers and T is the 
reciprocal of the sampling rate.) A low-price ageing clock oscillator can have a frequency offset of 50 ppm (In fig. A , 
this corresponds to a normalised clock frequency offset e=5.10- 5 ). Then at a Signal to Noise Ratio (SNR) of 20 dB, the 
25 ici degradation can represent more than 1 dB. Moreover, the parasitic rotation due to f can represent 1 5% of the angle 
separating neighbour points of the QPSK (Quadrature Phase Shift Keying) const ellatfon. The parasitic rotation due to 
e can represent 20% of this angle. Therefore, the error probability rises up severely. 

[0007] Known methods which exist for reducing clock and carrier frequency offsets work on the assumption that time 
synchronization Is achieved and that carrier frequency offset has been compensated up to a certain precision. However, 
so there Is still a need for a method with reduced complexity and better performance. 

[0008] In cable transmission , a sub-carrier Is dedicated to clock recovery, which Implies a lose of useful capacity of 
the system. 

[0009] In US_A~5345440 (Qledhill) a method is presented for estimating the carrier frequency offset and the clock 
offset. But only a formula for, straight carrier frequency offset estimation is provided, whereas^no formula for directly 

05 estimating the clock frequency offset Is given. Besides, carrier and clock are estimated separately. Bbth pbenomehona 
produce similar effects; so that a joint method for estimating both would have an Improved efficiency and a reduced 
complexity. .« l v • ..; .ft { . r 1 -\ .-.-u i .! 

, i [001 0] 1 1 n Elec. Lett 34(6) 520-521 March 1 998 , a method for jointly estimating carrier and clock frequency offsets 
is proposed. However,- this" method has the disadvantage of poor performance on a f reqUency-6 elective channel^ such 

40 as the radio-mobile channel or cable channels. j - » ■ - . ... • 

[0011] In US_A_5802117 (Ghosh), a joint method working on any channel is presented. However, It Is not a 'blind' 
technique le. It relies on aspecific training signal. Therefore, rt cannot be applied on any existing standardised com-: 
munication system. J 

[0012] There Is subsequently a need for jointly and blindly estimating caniBr and clock frequency offsets on any 
« channel, including frequency selective ones and rt is this need that the present invention addresses. 

[0013] Accordingly, the present invention comprises apparatus for estim atin g carrier frequency offset and sampling 
frequency off set between transmitter circuitry and receiver circuitry which communicate over a channel of an OFDrvT 
system, the apparatus including In a receiver circuit; 

so a focal oscillator having a frequency fb for converting a received OFDM modulated signal s(t). representing infor- 

mation symbols S(k) having components Sm(k) and sampled at a f irst sampling rate ft and modulated onto a carrier 
of frequency fc, to a base-band signal, 

an analogue to digital converter having a second sampling rate fr for sampling the baseband signal, 
a demodulator for performing a discrete Fourier transform on the sampled baseband signal to generate blocks of 
B5 symbols R(k) having components Rm(k) representing the information symbols S(k) of components Sm(k), and 

characterised by ; 

aflrst module lor removing parasitic effects of th e channel on R(k ) to g enerate at le a st one block Y( k) of components 
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a seco nd module for removing modulation effects from Y(k) to generate at least one block 2(k) of components Zm 
(k) and to compute V of components Vm where Vm^/Ym / , 

and a third module for estimating terms A and B rela ting respectivel y to the carrier f requency_oflset flc-fb/and the 
sampling frequency offset /ft-fr/ by performing a joint" maximun likelih ood ggtjrjrj^ " 

f * 
#»*— 

£[-V m y\(Z m )sin(A + mB) + K ffl 3(ZJcos(,4 + wfl)]= 0 
x 

K 

^X-mVJi{Z m )tm{A +mB) + mV m Z{Z m )cos(A + mfl)] = 0 

K 

-I 



[001 4] The Invention thus provides a method and appara tu s for Joint ly estimati ng carrier andc lock frequenc y offse ts 
on any ch annel, including Irequency selective one s. 

zo [0015] The invention computes the estimates f of /and e of,©, by analyzing the rotations produced by clock and 
carrier frequency offsets on the received signals. Tho estimate /can be added to that obtained with known algorithms, 
and may be used to compensate the carrier frequency offset on a time-domain signal. The estimate e can be used to 
steer the receiver clock frequency (e.g: using a Voltage Controlled Crystal Oscillator VCXO). Thus, the effect of carrier 
and clock frequency offsets is brought down to a negligtole level, at the expense of slightly higher complexity. 

25 [0016] Advantages of the method compared with existing techniques are : 

• The complexity required is low (a few times K, for a typical precision, where K is the number of useful sub-carriers). 

• The implementation of the invention allows either to reduce the cost of the receiver by relaxing constraints on the 
clock oscillator and on the carrier synchronization mechanism, or to save power (the gain can then be greater than 

• Since the method used is blind, It can bB applied on existing standards such as ADSL, DAB or DVB-T. It works 
either on coherent or on differential modulation schemes; 

• The method is flexible. The ability to track /and o depends on the number of symbols over which the estimation 
is performed. The more symbols, the more accurate and precise the estimate, but the lower the ability of tracking. 

3S 1 ; , " Furthermore} the symbols 1 picked up for the estimate do not need to be consecutlveiFor instance, In a time-varying 
channel; the chosen symbols can be spaced out. Thus, the complexltydecreases and the estimation takes benefit 
from diversity effects. To summarize, it is easy to make a trade-off between complexity, performance and tracking. 
5' •■; * ■ ■• +'X ><* ■ }, '-.-"■'-:*. ' .f • .=►";. »■.■■■ ■• • ,5 * . * ; 

; > [0017] V1 The lnvention can be used to advantage in the case of DAB because the method used is blind and can be 
*o easily applied on the differentially coded QPSK transmission; • . • r, • .h V; vV/s • '■ 

[0018] The Invention can-also be used as an arte rnative to the use of a dedicated, channel .>\ _ ,f * • 
- [0019] (Some embodiments of the Invention will now be described, by way, of example only, wfth reference to the 
drawings of which : - ; 

4s Fjg.r Is a graph showing number of sub-carriers v Signal-to-noise degradation due to ICi in an OFDM system 

having a relative clock frequency offset of 5.1 0exp(-5) ; 

Flg.2 Is a graph showing carrier Irequency offset v SlgnaMo-nolse degradation In an OFDM system ; 
Fig. 3 Is a schematic block diagram of a conventional OFDM transmitter; 

Fig.4 is a schematic block diagram of an OFDM receiver Incorporating carrier and clock frequency offset compon- 
50 sation means in accordance with the invention ; 

Fig.5 is a schematic block diagram illustrating the operation of the embodiment of Fig.4 in greater detail; 

Fig.6 is a graph showing sub-carrier Index v normalised frequency response for a typical OFDM channel; 

Flg.7 is a graph showing number of symbols v RMS clock frequency offset estimation error for the OFDM channel 

operating In accordance with the Invention ; 
gs Fig.8 is a graph showing number of symbols v RMS clock frequency offset estimation error for a DAB mode, OFDM 

channel operating in accordance with the Invention ; and 

Fig.9 is a graph showing number of symbols v RMS carrier Irequency offset estimation error for a DAB mode, 
OFDM channel operating in accordance with the invention. 
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[0020] With reference to Flg.3, an OFDM transmitter comprises an Inverse Fast Fourier Transform modulator 1 r a 
guard Interval Insertion block 2, a Parallel to Serial converter 3 a Digital to Analogue converter 4 and a mixer 5 having 
alocaJ oscillator input on line 6. Information symbols S(k) are fod to the Input of the modulator 1. The outputs from the 
modulator pass through block 2 where a guard Interval Is Inserted. After guard Interval Insertion, the signals are fed to 

5 the Parallel-to- Serial converter 3 and then the resulting multiplexed output is fed to the Digital to Analogue converter 
4. The Digital to Analogue converter 4 has a sampling frequency ft = 1/T which Is also the total bandwidth of the 
multiplexing process.The resulting analogue signal Is modulated onto a carrier of frequency fc by the action of the 
mixer 5 tor transmission over a radio channel to the receiver of Flg.4 There are N sub-carriers with K useful sub-carriers. 
The bandwidth is split among the N sub-carriers, so that the sub-carrier spacing Is-^. 

io [0021] Refering now to Fig.4, an OFDM receiver comprises a mixer 7 for converting the incoming signal on line 8 to 
baseband. A baseband signal at frequency 1b,( nominally the carrier frequency ) is supplied to the mixer 7 by a local 
oscillator 9. The baseband signal is then sampled In an Analogue to Digital converter 1 0 at a sampling rate f r (nominally 
the transmitter sampling rate, ft). 

[0022] The difference between the carrier frequency tc at the transmitter and the local osci llator frequency fb at the 
,5 receiver is called the carrier frequency offset (fc^fb). Expressed In terms of the number of sub-carrie r spacing s, *=NT 
(fc-fb) The diffe rence between the sampling frequenc y at the receive r and a t the transmitter is called the clock frequency 
offset (fr-ft) . The relative (or normalised) clock frequency o ffs et is denoted " " 



[0023] The output signal from the Digital to Analogue converter 10 is fed to a Serial to Parallel converter 11. The 
resulting de-multiplexed outputs have the guard interval removed by the guard interval remover block 1 2 and then pass 
to a Discrete Fourier Transform demodulator 13. The outputs of the Discrete Fourier Transform demodulator 1 3 are 
subsequently passed through, in sequence, a Remove Channel Module 14, a Remove Modulation Module 15 and a 
oo Computation of Estimates Module 1 6. The outputs of the Computation Module 1 6 are the desired estimates of /and o. 
[0024] The demodulator 13 produces blocks of symbols Rm(k) at Its outputs. In order to generate the estimates of 
/and e, the apparatus of Fig. 4 measures the phase drift due to carrier and clock frequency offsets between symbol 
(kr1) and symbol k. However, the phase has drifted continuously during the previous k-2, k-3. ... symbols because of 
. carrier and clock frequency offsets and also because of the effects of the channel.. The Remove Channel module 14 
. . 35 : j>k isolates th a phase drift du a to carrier, clock and chann e I between symbol (k-1 )and aymbo 1 1 k. « 
>= k •• [0025] :■■■! Hence, at the output of the Discrete Fourier Transform (DFT) * demodulator 1 3, the A** block of symbols 
*i ; V k% ^(k) - % m.£$,m * 0 Is processed by the Remove Channel Modulo 14. Here, the parasitic effects of the 
> channel and of the rotation due to carrier and clock offset on block k - 1 are removed. At the output of. the Module 14, 
j. : . tr ; i the block of symbols Y^fk) -{f S m £-£.m * 0 is obtained. The processing procedure^ Module 14 can consist either 
- 4o in differential demodulation or In equalization. For instance, in DAB, we have Y m (k) = tf m (k)Ffo(k-1), > 



[0026] A second processing step In the Remove Modulation Module 1 5 removes the "useful* rotation , due to mod- ». 
v ulatlon, from Y M (k) and also computes the quantity V n ;. a i vy. Thus there, remains only the parasitic rotation produced 
by carrier and frequency offsets. In addition with various sources of perturbation (including thermal noise). The output 
of the Module 15 is the block Z m (k) -fsm S$,m * 0. For example, in DAB, the rule can be the following: 



20 




2S 




45 



• if* (VmW) > 0 and 0 (V^A)) > 0 then 



50 



• it 9t { Y m {Kj) < 0 and 3 ( Y^A)) > 0 then 
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« * [Yjm < 0 and 3 [YJK» < 0 then 
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• « 91 ( V m (*» > 0 and O ( K m (*)) < 0 then 
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But in the case of DAB, the rule can also be: 12^(^1 = \Y m {fy\ argr(2^(A)) = 4*a/g( V m (A)) 

[0027] In a third step, a joint Maximum Likelihood Estimator (MLE) tor (r» la computed from the observation Z m (A) 
-£ <. m s£,m * 0. The Joint MLE requires the resolution ol the following system of equations: 




1S 




K„JR(Z w )sin(,4 + /wfl) + K M 3(Z.)cos(^ + «J?)]= 0 
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)sin(/4 + m*) + mK w 3(Z w )cosM + mB))= 0 

2 




2$ 






30 




(equation 1) 






It can be verified that, as long asl4l + 




■ i-U- :■ ... 1; . 

• ■»■..■ ,c ■' 3$ - >: » 










is small compared toil , a I 8 ' order approximation of the sine and cosine terms is enough. Then/the third processing . 
step consist* in cbmputlng*the following 'coefficients :y * ... > ■ •. i ..■ 








,^ , .... ..... . . ..... ....... . ... .. 
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«, = 2«K„SR(ZJ ; a, = J t m i V m 'W m ) > «« - £^91(Z„,) 
_f ^ 

*r*0 **0 
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= £rnV m 3{Z m ) ; a, = ^.3(2.) 
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a>a«-aj 

[0028] A and 0 are proportional respectively to the desired f e and & The ratio depends on which method has been 
10 chosen tor removing the modulation. For example, in DAB, if the first rule was applied, then : 



15 



f 



N 



E — ■ 



2n(N + D) 
2n(N+D) 



20 [0029] But if the second rule was appli ed, then: 



o^TV + D) 

30 Where N Is the number of samples and D is the length of the guard interval in terms of sample number 
The MLE system of equations Is derived as follows. 

By the action of the Remove Channel Module 1 4 and the Remove Modulation Module 1 5, the vector 

.. . is obtained.. - . . . - . .-. .■ 4,3 v-- - , \ 

This vector Is approximated to the model: =' : r , ■ ' : 

40 . * • ■ . . ■■ ■ . v ........ . _ . ' . 

; ^ *m<*) = VJkjgJ"? + 7 m (/f) With - £ 5 m.S.£>*0 

where 

45 ^m^O represents the modulus of the useful part of the signal, 
V m (A) = \Z m </r)l = I >y*}l for high eignai to noise ratio.. 
y:AWQN (Additive white Gaussian noise) 
lit Is desired to fit /and ©so that the model best fits reality. 
Introducing vector 

so 

*2 2 

ss the Maximum Likelihood Estimator 

tor f and e according to the observation Z is the mathematical formula: 
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(fej)= ArgMax iJc) [-(Z-Ve»e jmB ) H (Z - Ve^e'"*)] 

s [0030] From which equation (1) above is derived. 

[0031] In an alternative embodiment, the estimation may be improved by using Lblocks 2(/c,), Z(k^ t ..„ Z{k$ instead 
of Just one block. 

[0032] In this case, the coefficients to be computed for the estimate are the 

10 

«f ,ie{1, 2,3,4,5} 



[0033] The total number of blocks chosen would depend upon the processing time available. 

15 [0034] The estimation may be performed on successive symbols kj = k1, k2 kL, with k2=k1+1, k3=k2+1 etc. 

However, when the channel variations are slow, the channel attenuation is almost the same on consecutive symbols. 
So, if channel fading worsens the estimate, using successive symbols will increase the complexity but not significantly 
Improve the estimate. In this case then, the symbols chosen for the estimation can be spaced out, for instance choosing 
k2=k1+5, k3=k2+5...etc. The spacing chosen will depend on how fast the channel is varying. 

20 [0035] tn a further alternative embodiment and in order to reduce the complexity of the process, coefficients can 
be replaced by 1. No loss of efficiency has been experienced in simulations. 

[0036] Figure 5 represents the processing architecture for realising the corresponding third processing step in which 
Vm is set to 1 (unity). 

[0037] In any of the embodiments, if the 1 st order approximation of cosine and sine is not sufficient, a second or 
& higher order approximation can be applied. 

[0038] The system of equations obtained then becomes non-linear and iterative resolution methods can be applied. 

[0039] Figs. 6 to 9 demonstrate the performance of the invention lor the cases of two types of communication channel. 

The normalised channel frequency response in the first case is shown in Fig.6. This corresponds to a baseband gable 

system. The second case relates to a typical DAB urban channel. For cases, N=512 sub-carriers, K=384 useful sub- 
30 carriers. For each case shown, c is random uniform In [- S.IO-^S.IO- 6 ]. in the first case, f- 0. In the second case, r*is 

random uniform in [- 3.10*3. 10^|. 50 simulations are performed for each SNR In {BdB,16dB;24d8}. The Root Mean 

Square (RMS) of the estimation error is computed and plotted. 

[0040] In the first case, the RMS of the clock frequency offset estimation error Is shown in fig. 7. Averaging only on 
2 successive blocks is enough to have a clock frequency offset below "lO.ppm. ...->. 
3S [0041] In the second case, the simulated system corresponds to DAB mode 2 on a typical urban channel. Here, the i 
*■ estimation is performedon 1 symbol out of 3, to take benefit from diversity effects. After 10 iterations, the clock frequency , 
« - offset Is below 20 ppm. The RMS of the clock frequency offset estimation error Is shown In fig. S. and the RMS of the.: 
carrier frequency offset estimation error in fig. 9. Even, in such poor conditions of transmission (fast, varying channel), 
- the method swiftly reaches the targeted performance . .. — : ; 

40 - 1 

' • r \ 

■' ^ ' Claims 1 - .* ." ■■}■■• .. V ; v .. - ; .. - 

1. App aratus for estimating carrier frequency offset and sam pling frequency offset betwe en transmitter circuit ry(1 S) 
45 - and receiver circuitry (B-16) which communicate over a channel of an OFDM s ystem, the apparatusjnduding jn 
a receiver circu it; * • " 

a l ocal osclDator (9) having at requency fb for converting a received OFDM m od ulated signal s(t), representing 
Infor matio n symbols S(k) having components Sm(k) sam pled at a first sampling rate ft arid moduiatWomo a" 
60 carrier of frequency fc, to a base-band 6i gnal> 

an analogue to digital converter (10) having a second samplin g rate fr tor samp ling th e baseband signal , 

a de modulator (13) for performing a discre te Fourier transform on the sampled baseband signal to generate 
'bjocks of s ymbols R(k) having components 1%i(k) representing the information symbols" S(k) of components 

Sm(k), and chara cterised by ; 

55 a first module (14) for removing parasitic effects o f the channel on R(k) to gene rate at least one block Y(k) of 

component s Ym (k). 

a second module (15) for removing modula t ion effects from Y(k) to generate a t least one block Z(k) of co m- 
ponents Zm(k) and to compute V of components Vm where Vm=/Ym /, 
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25 



30 



40 



45 



50 



a nd a third module (15) for estima tin g terms A a nd B relating respectively to the carrier frequency offset he- 
tb/_ond the sampling frequenc y oflael/ft-1r/ by performing a Jo int maxlmun likelihood estimation (MLE) such that; 



10 



15 



£ [- VJf{{Z m ) sin04 + mB) + V„3(Z m ) cos(^ + mB)] = 0 



<n*0 
K 



2! [- m yjto&m ) + rnB) + m V m 3(Z OT ) cos(A + mB)] m 0 



20 



2. Apparatus according to claim 1 in which the second module (15) generates a plurality of blocks Z(k) and the third 
module (16) adjusts coefficients of the MLE accordingly. 

3. Apparatus according to claim 1 or claim 2 in which Vm is set to unity. 



. 35 
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